Theoretical study of X-ray absorption of three-dimensional topological insulator Bi 2 Se3 
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X-ray absorption edge singularity which is usually relevant for metals is studied for the prototype 
topological insulator Bi2Se3. The generalized integral equation of Nozieres and Dominicis type for 
X-ray edge singularity is derived and solved. The spin texture of surfaces states causes a component 
of singularity dependent on the helicity of the spin texture. It also yields another component for 
which the singularity from excitonic processes is absent. 
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Introduction- The topological insulators (TI) have been 
intensively studied recently^ 2 - TI have bulk energy gap 
but they have conducting (namely gapless) states at 
boundary^ Quantum Hall state is an example of 2- 
dimensional TI with broken time reversal symmetry, and 
the conducting states at boundary are nothing but the 
well-known edge states^ There exist 3-dimensional TI 
with time reversal invariance, and they have conducting 
surface states (SS) which are protected by Z 2 topolog- 
ical invariants in the bulk4r— The energy band of the 
SS takes the form of Dirac conei^ SS have been first 
observed in Bii-^Sba;, but many of their important fea- 
tures were not clearly discerned due to small bulk gap and 
disorder effect. 7 Stoichiometric TI possessing the sim- 
plest SS structure, namely a single Dirac cone, have been 
proposed for Bi 2 Se3,Bi 2 Te3,Sb2Te34i^ The single Dirac 
cone SS has been observed in ARPES experiment for 
Bi2Se3fiS These materials can be realized as TI owing to 
the band inversion mechanism driven by large spin-orbit 
coupling^. The spin texture which is a distinguished fea- 
ture of SS has also been observed in the spin-resolved 
ARPES experiment fii 

X-ray absorption (and emission) spectroscopy is a very 
important method in the study of the electronic structure 
of core electrons. An incident x-ray photon excites a deep 
core electron to an unoccupied state with higher energy, 
leaving behind a positively charged core hole which can 
be treated to be immobile in many casesi 12 i 13 If there 
exist conduction electrons (of metals), they react to this 
suddenly created potential by deep core hole. The con- 
duction electrons interact with deep core hole in two dis- 
tinctive ways: the excitonic process 14 which is essentially 
attraction between hole and conduction electrons and 
the orthogonality catastrophe 15 which means a vanishing 
overlap between the ground-state wavefunctions before 
and after the creation of the deep core hole. Both exci- 
tonic process and orthogonality catastrophe are singular 
near x-ray absorption edge for fermi liquids, and they re- 
quire nonperturbative treatments For TI, partially 
filled SS comprise the conduction electrons in spite of 
energy gap in the bulk. Evidently, it is of interest to 
investigate how the above singular behaviors for the con- 
ventional fermi liquids are modified for the conducting 
states realized by SS of TI. 

However, there is a caveat. In the X-ray absorption 



experiment with photon incident perpendicular to sur- 
face, X-rays penetrate deep into the bulk of a sample, 
hence providing its bulk properties^ The SS of TI re- 
side near surface, so that in this experimental setup the 
contribution to the absorption from SS is expected to be 
rather small. The decay length (along surface normal) 
of SS of Bi2Se3 can be estimated to be in the range of 
4 - 10 A using Eq.(32) and Table IV of Kef. !) . For 
the substantial amount of X-ray absorption to take place 
in conjunction with SS the attenuation length of X-ray 
should be comparable to the decay length of SS. The at- 
tenuation length of X-ray can be controlled by its energy 
and the incident angle (measured from surface). Taking 
various factors mentioned above into account, we choose 
to focus on a core level N 3 4p 3 / 2 of Bi whose binding en- 
ergy is 678.8 eV. At this energy the critical angle is 3.26 
degrees and the attenuation length at 3.4 degrees is 40A 
which is indeed comparable to the decay length of SS. 
We note that the attenuation length for the normal inci- 
dence is about 1,000 A, 2 2, Presently there seems to be no 
experimental report on the X-ray absorption in Bi2Se 3 . 

In this Brief Report we report the results on the X-ray 
edge problem of Bi2Se3. The spin texture structure of SS 
modifies the singular edge behavior compared to that of 
conventional Fermi liquids. The most salient differences 
are the appearance of a contribution depending on the 
helicity of spin texture and the other which is free of the 
singularity from excitonic process. The main result of 
this Brief Report is Eq. (f34|) . 

Setup- The Hamiltonian for the SS is given by (see 
Eq.(34) of Ref.0, and a, /3 =f, J, denote spin ) 

Hss = ^2 4 Q h Q/3 c k(3 , h = (h Q/3 ) 

a,/3=t,4- 

h = (C + C 2 (k 2 x + k 2 y ))l 2 + A(a x k y - a y k x ), (1) 

where a x>y are Pauli matrices acting on spin space (I 2 
is a unit matrix) and Ck Q is the destruction operator 
for SS with wavenumber k and spin a. In Eq.([T]) we 
have ignored the trigonal distortion terms proportional 
to {k x ± iky) 3 . The numerical values of the parameters 
of Eq. dTJ which are appropriate for Bi2Se3 are given by 

Co = 3.37x 10- 2 eV, C 2 = 23.7eV-A 2 , |i| = 3.30eV-A. 
The Fig. 4 of Ref.[9[ suggests the wavenumber cutoff be 



2 



0.1 A . The energy eigenvalue is given by 



E ± (k)=Co + C a fcfi±|A|fc,|, h = Jkl 



(2) 



When the C 2 fc 2 is much smaller than \A\k\\ the Dirac 
cone structure with apex at k = is manifest. 

The non-interacting Matsubara Green's function of SS 
can be expressed as (the hat denotes matrix and <f> is the 
azimuth angle in x-y plane) 

g(ie, k) = l 2 96 + [(+i)Si2e-** + (-tj^ie**] g a , (3) 

where is a 2x2 matrix whose only non-vanishing ele- 
ment is 1 at entry, and 



fld(ie,k) = 



ie + fi — Cq — C* 2 fc 2 



[ie + n - E+ (k)] [ie + fj,-E_ (k)] ' 



3o(«e,k) 



[ie + /x - (k)] [ie + n-E- (k)] ' 



(4) 



where /i is the chemical potential, and fj, > Co will be 
assumed (namely, the lower Dirac cone is completely oc- 
cupied). The interrelation between the spin and the an- 
gle in Eq. © is nothing but the manifestation of the spin 
texture of SS. 9 Note that the Green's functions of Eq.fU 
are independent of the angle 4>. The Green's functions 
summed over wavenumber in the long time limit are given 
by (r is imaginary time) 



9a{t) = , 9o(r) 

T 



sgn(l) ff d (r), 



(5) 



where p = Akp/ Attvf is the density of states at Fermi 
energy (A is the area of unit cell), and vf and hp is 
the Fermi velocity and Fermi momentum, respectively, 
whose detailed form does not concern us here. Note the 
sign factor sgn(A) in Eq.©, which is the signature of the 
helicity of the spin texture of SS.— 

The core level N 3 4p 3 y 2 of Bi is labeled by the z- 
component of the total angular momentum J = 3/2. 



H 



h..le — Eh 

mj=±3/2,±l/2 



ft* b 

u mj u mj j 



(6) 



, is the creation operator of the hole. Eh is 



where bj, 

the core level energy, and \x + Eh is the (unrenormalized) 
threshold energy for X-ray absorption. The potential cre- 
ated by the deep core hole will be assumed to be spher- 
ically symmetric, and for simplicity we will consider the 
isotropic scattering only, so that the potential scattering 
matrix clement for SS is simplified tc 18 



-Vb, Vq > is constant. 



(7) 



Then the interaction Hamiltonian between SS and the 
deep core hole is given by 

Hint = E(- y °)(£ <Lc*«)(E h lM,i (8) 

k,k' a mj 



where a suitable cutoff in the wavenumber sum is as- 
sumed implicitly. The total Hamiltonian consists of 



H, 



H 



ss 



H 



hole 



(9) 



With the Hamiltonian Eq.© the hole quantum number 
mj is conserved, and it implies that the deep core hole 
Green's function is given by 

- (brnA^bly)) = 5 mj ,m>p{T - r'), (10) 

where the function D{t) is independent of mj. The X- 
ray absorption intensity I(lo) {lu is the frequency of the 
incident X-ray) can be expressed in terms of correlation 
function [using Eq. (fTU|) ] as follows: 21 



I(u>) = Im 



F ( r ) = E E M ^( k ' a i m ^) M qA(k',^|my) 

k,k',a,/9 mj 



x F, 



Fw a p\ mj (T) = (cka(r)& mj (r)6t lJ (0)4^(0)), (11) 

where M q ^(k, a\mj) is the X-ray transition matrix el- 
ement from the deep core state \mj) to the SS (Bloch 
state) |k, a), and q and A is the wavenumber and the 
polarization of the incident X-ray, respectively. In many 
cases of interest, the wavenumber (q, k) dependence of 
the X-ray transition matrix element can be ignored. This 
is due to the localized nature of the wave function of core 
electron. In the presence of the strong spin-orbit cou- 
pling such as the case of Bi, the electron-photon interac- 
tion receives additional contribution from the spin-orbit 
coupling. 22 The explicit form of the transition matrix el- 
ement including spin-orbit contribution is 

M qA (ka|mj) = J dVe^l^^-^eqA-W™.,] 

/M 6 Wk l? • [ - ^e qA X V* mj + W X e q A*m. 7 ] }, 
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(12) 



where e q A is the polarization vector of X-ray and V(r) 
is the periodic crystal potential. 4>k{r) and ^> mj (r) is 
the (spinor) Bloch wavefunction of SS and the (spinor) 
wavefunction of core electron, respectively. The dipole 
approximation e lq r ~ 1 will be assumed below. 

Correlation functions- The correlation function 
Fkk'a/3\mj( T ) of Eq.(JTT} can be obtained from 

^kk' 

Tl,T%) — (TcUa{0^rn J (Ti)b mj (T2)c k , l3 

(£')> 

(13) 

by the limiting procedure £ — > t\ — t c and £' — > t 2 + t c (t c 
is a short-time cutoff). T denotes time ordering. For the 
absorption we have to take T\ > r 2 . We employ the equa- 
tion of motion method of Ref.pjl to derive the integral 
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equation for Fkk'Q/3|m 7 (^C'l r i: T 2)- The conservation of 
the hole number causes the equation of motion to close 
on itself,— and we find the integral equation 

Fw a /3\ m j(£,t'\Ti,T 2 ) = 6i<k>g a p(€-£\k)D(Ti -r 2 ) 

rri 

+ X] / dT 9 al (£ - T, k)Vkq'-F q / 

k'7/3|m j 

(t, £'\ti,T 2 ), 

i J To 

q',7 2 

(14) 

where g a p(£, k) is the Green's function Eq.Q in (imagi- 
nary) time domain. Decomposing F]^i a m mj as follows 

FuU> a/3\mj(£,A' \ T 1, T2) = Gkk'a/3 (£, £' I T l , T 2 )D{t\ - T 2 ) 

(15) 

Eq. (fl4l becomes the following integral equation: 

drg a7 (^ - t, k)y k q'G q / k ' T ^(r, £'\n, r 2 )(16) 

q-, 7 2 

Eq. dTSl) is the generalization of Eq.( 17a) in Ref.[18] to our 
case of SS. G k k'a^ and D{t\ — r 2 ) of Eq.tfTB")) represents 
the excitonic processes and the orthogonality catastro- 
phe, respectively^ It can be shown that the hole Green's 
function D(t± — t%) can be obtained from the solution of 
Eq. ([T^l) via parametric integral (see Eq.(21) of Ref.[l8j 
and Eq.(ll) of Ref.Ql]). Thus once Eq.® is solved, we 
can find X-ray absorption intensity from Eqs. (jllll5|) . 

In fact, we need to find the Green's function Gkk'a/3 
summed over wavenumber weighted by transition matrix 
element M q >(ka|mj). In most cases of simple metals the 
wavenumber dependence of transition matrix element is 
ignored. However, in our case such dependence is cru- 
cial because, as can be seen in Eq.©, the spin texture 
structure is encoded in the angle dependence of Green's 
function. In view of this we expand the transition ma- 
trix element in Fourier series of e 1 ^ but we will ignore 

fc|| = y/fc^ + k% dependence. 

M qA (ka|mj)« J2 e in *M^(a\mj), (17) 

71 = 0, ±1 

where only n = 0, ±1 terms are kept since the higher 
order contributions will be smaller because they involve 
higher power of fciir. Let us define (henceforth time ar- 
guments are suppressed for notational clarity) 

G a p\n,n' = X] e m( ^(e m * )*Gkk'a,3, (18) 
k,k' 

where a cutoff in wavenumber sum is implicitly assumed. 
Applying the definition Eq. flTS)) to Eq.© we find 

9a0\O,O = 9ap\l,l = 9a0\-l,-l = <W <7d CO > 

9a/3\i,o = 9 a p\o,-i = S a ^dpi.(+i)g (T), 

9a0\O,l = 9a0\-l,O = 8ai.80-t(—i)g a (T), 
9a0\l,-l = 9a0\-l,l = 0. (19) 



Now Eq. (fl"6]) can be recast into the following form [ recall 
Eq.0 ] 

G a 0\n,n'(£> C\ T li T l) — g a 0\n,n'(£ ~ O 

+ (-Vq) [ dTg ajlnfi (^ - t)G 7( 3| ,„'(t,£'|ti> t 2)(20) 



which are coupled integral equations. Noting the factor 
9a-y\n,o an d from Eq. (jl9[) . we find the nontrivial solutions 
(with non- vanishing 2nd term) obtain only for n — - 0, ±1. 

Solution of integral equation- The well-known Nozieres 
and Dominicis (ND) (asymptotic) solution in long time 
limit is [g(£) — — |=non- interacting Green's function] 18 



G N D(^'|n,T 2 ) = cos 2 5g{€-?) 



(21) 



.(n -£)(£' 

which satisfies [ compare with Eq. ([2H|) ] 
G N D«,e'|r 1 ,r a )=^-0 
+ (-V ) r drgfc-^GimfaZ'ln,^), (22) 



where S is the s-wave scattering phase shift 

6 = tajr 1 [irV p]. (23) 

For n = n' = 0, using Eq.([19]), Eq.® is found to 
reduce to Eq. (|22|) . Hence (time arguments suppressed) 

G a( g|o,o = <WGnd- (24) 
Noting Eqs. (|5ll9p we also find the solutions for n' = ±1 

G Q/3 |o,i = 5 ai ,5^(-i)sgn(A)G N T>, 
G a 0\o-i = 6 at S/3i(+i)sgn(A)G ND . (25) 

For other values of n', G a p\ 0n i = 0. 

Next consider the case of n = 1. From the property of 
SWyli.o [ see Eq.(|19p] the nontrivial solutions obtain only 
for a =t, so that 



G\.p\\,n' — 9\p\x,n>- 

Now take a =\. 

G T/ 3|l,„'(£;£'l r l; r 2) = 9f0\l,n'(^ ~ O 



(26) 



+ (-Vb) / drg nll ^-T)G imnl (r^'\ri,r 2 ).(27) 

If n' = —1, then from g~fp\i,-i = and Gj^iq,-! = [see 
Eq.([25|)]. we conclude that 



G 



t/3|l,-l 



0. 



(28) 



For the case of n' = of Eq. (j2"7| . (3 should be I, otherwise 
both g and G vanish. Thus 



G 



ttlio 







(29) 
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and for j3 =\., multiplying both sides of Eq. (l27j) by 
(— i)sgn(A), we find the equation becomes exactly ND 
type Eq.(|%2]l. so that 



G n \ w = (+ l )sgn(i)G ND . (30) 
Repeating the similar analysis for other cases we obtain 

G Q /3|0,0 = da/sG^D, G Q ^| 1: _! = G Q/ 3|_ 1:1 = 0, 



r , _ i Gnd , , _ 

L-a/3|ll - \ q ], Cr Q/3 |_i_i - 



3d 
Gnd 



G a /3\io — G Q ^|o-i 



0\ 

-zsgn(i)G ND Oj ' 

isgn(A)G N D 




(31) 



The leading behavior of the core hole Green's function 
D(t) can be obtained from the second order linked- 
cluster expansion^ The important contribution turns 
out to be 



d( / d£'Tr 



5 (o) (e _ n5 (o) ( ^ 



, (32) 



where gW (£-£') = £k s( k > £ - O = <7d(£ - £%• Thus 
the helicity of the spin texture of SS does not affect the 
core hole Green's function, so that the situation becomes 
essentially identical with that of ND solution. The evalu- 
ations of the integral of Eq. (|3"2")l yields logarithms, which 
is to be exponentiated in the linked-cluster expansion. 
Comparing with ND solution we find (N c = 2, (t, 4-)) 



D{t > 0) - e 



1 



(r/r c )^(V-) 2 



(33) 



is the renormalized threshold for X-ray absorption. 
Results- Combining the solutions Eqs. (|31l33[) with the 
transition matrix element Eq. (|12[) we can obtain the re- 
sult for X-ray absorption intensity. Eq. (I12[) has two com- 
ponents: one from direct dipole transition and the other 
from spin-orbit coupling. Dipole transition conserves the 
spin, so that it is diagonal in spin (namely, M a M^ cx 5 a p. 
This can be verified explicitly in our case). Then Eq. (|3"Tj) 
tells us that only G Q/ 3|(oo),(ii),(-i.-i) contribute. Among 
these, Gnd includes the singularity from excitonic pro- 
cesses while gd does not. The cross term of dipole tran- 
sition and the spin-orbit contribution allows spin off- 
diagonal configuration. From Eq. (|3"Tj) we find that this 



contribution is proportional to the helicity sgn(^4). The 
above considerations give (Q(x) is step function) 



+ (d + sgn(i) Co )( W - lo* t )- 2S ^+ n ^ s I^ 2 ],{M) 

where c d ,Cd,c are constants. 

Summary and concluding rcmarks-We have studied the 
X-ray absorption edge singularity of the prototype TI, 
Bi2Se3. For the singularity to exist, gapless conducting 
states are necessary, and SS provide those. Due to the 
spin texture of SS, two interesting modifications com- 
pared to that of convention metals arise: (1) helicity de- 
pendent contribution [sgn(A) term of Eq. (|3"4")l ] (2) con- 
tribution free of singularity from excitonic process [ the 
first term of Eq. (|34|) ]. These features can be verified ex- 
perimentally by extracting the surface contributions us- 
ing the angle and energy dependence of the penetration 
depth of incident X-ray, and we have also suggested a 
specific core level appropriate for experiments. The an- 
gle and energy dependence of the penetration depth can 
be also used in distinguishing the helicity effect of TI 
from that of the conventional Rashba spin-orbit energy 
bands since the latter reside in the bulk. We mention 
in comparison that the graphenes have two Dirac cones 
from valley structure, so that its qualitative properties 
different from those of TIi^ 3 . We also mention that when 
the Fermi energy crosses the apex of Dirac cone, all of 
the singularities disappear. 23 * However, this situation is 
not generic for TI and is not elaborated in this Report. 
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